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Abstract 

We study a set of textures giving rise to correct masses and mixings of the 
charged fermions in the context of leptogenesis. The Dirac neutrino texture 
pattern is assumed to be identical with the up quark texture. The heavy 
Majonana neutrino mass matrix is obtained by inverting the type-I see-saw 
formula and using the neutrino masses and mixings required by the solar and 
atmospheric neutrino oscillation experiments as input. After making the fea- 
sibility study of the generated lepton asymmetry via the decay of the heavy 
right handed neutrino, we compute the generated baryon asymmetry by nu- 
merically solving the supersymmetric Boltzmann equations. We find for these 
models that both the hierarchy of the texture as well as the placement of the 
texture zeros are important to the viability of leptogenesis as the source of 
the observed baryon asymmetry of the universe. 
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I. INTRODUCTION 



Relationships between fermion masses and mixings have been the subjects of much the- 
oretical interest starting with a postulated relationship between the Cabibbo angle and the 
down and strange quark masses. Most of the unknown parameters in the Standard Model 
(SM) occur in the Yukawa sector and any relationships between these parameters is welcome 
theoretically and experimentally testable. The interest in models of fermion masses and 
mixings accelerated with the advent of grand unified theories. In these models, the gauge 
multiplets of the Standard Model are unified into multiplets of the grand unified gauge group, 
and relationships between the parameters emerge naturally as a consequence of the larger 
symmetry. These models can be augmented by global symmetries or texture zeros in the 
Yukawa coupling matrices can be assumed, leading to further predictions. 

A mystery of the Yukawa sector is the obvious hierarchy that exists in fermion masses 
and mixings. The top quark mass is much larger than the charm quark mass which is still 
much larger than the up quark mass, for example. The (Cabbibo-Kobayashi-Maskawa) CKM 
matrix is measured experimentally have small mixing angles. Clearly a fundamental theory 
that explains the origin of the couplings in the Yukawa sector rather than just parameter- 
izing them, should explain these features. The Yukawa sector of the Standard Model is 
parametrized in terms of 3 x 3 matrices, so the hierarchy exhibits itself as a hierarchy among 
the elements of these matrices. 

The usual predictions from these models of fermion masses and mixings are relationships 
between the masses of quarks and leptons or between the mixing angles of the CKM matrix 
and dimensionless ratios of the quark and lepton masses. In addition there are often pre- 
dictions for the amount of charge conjugation-parity (CP) violation and predictions for the 
CP asymmetries of meson decays. It is of interest to consider whether further constraints 
are obtained after making reasonable assumptions about some other physical observable. In 
this paper we consider textures (or patterns of zero entries) of Yukawa coupling matrices and 
assume that the baryon asymmetry has its origins in the decay of heavy Majorana neutrinos 
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which violate lepton number. The asymmetry in lepton number is recycled into a baryon 
number asymmetry via the sphaleron process. This idea was first put forward by Fukugita 
and Yanagida |IJ and has come to be known as baryogenesis via leptogenesis. While this 
constraint is admittedly more speculative than the comparison of masses and mixing an- 
gles derived from experiment, it is instructive to determine which properties of the Yukawa 
textures are essential for the baryogenesis via leptogenesis to work. 

Sakharov pointed out || that a small baryon asymmetry may have been produced in 
the early universe if three conditions are satisfied: 1) baryon number is violated, 2) charge 
conjugation symmetry (C) and CP are violated, and 3) there is a departure from thermal 
equilibrium^. Since both C and CP are violated in the Standard Model, and baryon number 
is violated by a nonperturbative effect called sphalerons, the natural place to look first to 
explain the generation of a baryon asymmetry is within the Standard Model itself. However 
this line reasoning does not work: the required Higgs masses is too small and has been ruled 
out by the direct searches at LEP [|J. One can try to extend the Standard Model: one 
popular attack is to consider the Minimal Supersymmetric Standard Model (MSSM) and 
assume that the source of the CP violation is still contained within the CKM matrix ||. 
One can achieve the observed baryon asymmetry, but only at the expense of going to a corner 
of parameter space, namely one requires a light scalar top quark (top squark or stop). This 
kind of solution should rightly be regarded as unnatural on theoretical grounds, although 
it is of much experimental and phenomenological interest primarily because it is just out of 
reach. 

We pursue in this paper the alternative approach of baryogenesis via leptogenesis: With 
the confirmation of the atmospheric neutrino anomaly in the SuperKamiokonde |5[], it seems 
plausible to start from the point of view that it is likely that the neutrino mixing is occuring 



1 A baryon asymmetry could arise even in thermal equilibrium if CPT is violated. See, for example, 
Ref. §. 
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and the neutrinos have masses given by a see-saw mechanism. The heavy Majorana neutrino 
mass matrix is obtained by inverting the type-I see-saw formula 

m cS = m T D ^—m D , (1) 
M N 

and using the neutrino masses and mixings required by the solar and atmospheric neutrino 
oscillation experiments as inputf]. The neutrino sector contains a new source of CP violation; 
the interference between tree-level and one-loop contributions to the Majorana neutrino 
decays can give rise to a lepton asymmetry. In this scenario, the amount of CP violation that 
gives rise to a lepton asymmetry and ultimately to a baryon asymmetry depends critically 
on the Dirac mass matrix of the neutrinos in two ways: (1) the overall hierarchy pattern of 
the matrix and (2) the placement of the texture zerosQ. Furthermore, the generated lepton 
asymmetry can be erased by subsequent lepton-number violating scattering and this dilution 
can depend on the placement of the texture zeros. 

In this paper we start from the position that the positive observation in the solar and 
atmospheric neutrino experiments suggest that there is a new scale of heavy physics. We 
assume that heavy right-handed neutrinos exist and the lightness of the observed neutrinos is 
the result of a seesaw mechanism. In this framework we study the plausibility of leptogenesis 
in the case of neutrino Dirac mass matrices with texture zeros and hierarchical structure 
similar to the ones that are consistent with low-energy data in the quark sector. 

II. YUKAWA TEXTURES 



2 We neglect the possibility that there are contributions from a left-handed triplet Higgs 

3 If a texture zero predicts a small level of CP violation at the GUT scale, this suppression will be 
preserved by the renormalization group scaling. 
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TABLE I. The list of textures suggested by RRR. 
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Ramond, Ross and Roberts (RRR) performed M a systematic search for all possible 
symmetric quark and lepton mass matrices with five texture zeros at the unification scale 
that are compatible with low-energy measurements. They found a total of five possible 
solutions, which we display again in Table H for convenience. 

We assume the Dirac neutrino mass matrix at the GUT scale has the same texture zeros 
as the up quark matrix 

m' D ~ m u — X u v sin (3. (2) 

In certain situations where the Yukawa interactions are minimal, grand unified symmetry 
enforces an exact equality. More generally one might expect the equality of elements in the 
neutrino texture and the up quark texture not be exact, but be related by Clebsch coefficients 
(very often 3). These Clebsch factors are typically small and do not upset the hierarchy of 
the matrices. The general qualitative features of our analysis is not affected by these factors 
of order one, since the amount of baryon asymmetry generated in a model with a particular 
texture is governed by the hierarchy (given in terms of the parameter A which is fixed by 
the Cabibbo angle) and the position of the texture zeros. 

III. NUMERICAL SOLUTIONS 

We consider small angle MSW solution of the solar neutrino problem through the mixing 
v e *-* v y, and maximal mixing solution of atmospheric neutrino oscillation through the mixing 

«-> v T . We take as inputs the following neutrino masses consistent with the experimental 
measurements, namely [§,0 

0.8 < sin 2 29 23 < 1 , KT 3 < Am 2 3 < 1(T 2 , (3) 
3 x KT 3 < sin 2 29 12 < 2 x KT 2 , 5 x KT 6 < Am 2 2 < KT 5 . (4) 

The inverse neutrino mass matrix is 
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/1/mi \ 
m^ 1 = l/m 2 

v o o 1/7713 y 

This can be rotated by a mixing matrix 

K = V 13 V 12 V 23 , 



(5) 



(6) 



where Vij is a rotation matrix in between the i,j rows and columns by an angle For 
example 

/ C12 S12 0\ 

-sia C12 

V o 01/ 

where Sy = sin 9y and Cy = cos 6ij. We have taken the mixing matrix to be real for simplicity. 
Then the Majorana mass matrix in this basis is 
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(7) 



M N = m'nVm-dV 1 m 



1t/T / 



D * '"•eff 



b D i 



The lepton asymmetry is created by the decay of the lightest of the heavy majorana neutrinos. 
Consequantly we have to go to a basis in which the Majorana mass matrix is diagonal. This 
can be diagonalized by a matrix K so that 



M^ g = K T M N K 



(9) 



Note the Dirac and Majorana mass matrices are related through Eqn ([!]). It can be easily 
seen that the Dirac neutrino mass matrix in such a basis is, 



mo — Km D K 



(10) 



The CP-asymmetries in the neutrino decays arize from the interference between the tree 
level and one-loop level decay channels 

1 1 



lnv 2 2 (mime).-,- n 



[w} D m D ) 2 nj 
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(11) 



where 



g(x) = y/x 



m 1 1±£'| + 



X — 1_ 

and v% — v sin /3. The other parameter of most interest is the mass parameter 



(12) 



rh x = {m °2 D)l1 , (13) 
Mi v ; 

which largely controls the amount of dilution caused by the lepton number violating 
scattering Q. A large enough lepton asymmetry can result only if fh\ is in the range 
10" 5 < fh\ < 10~ 2 [13|. For too small values of rhi, the Yukawa interactions are too 
weak to bring the neutrinos into equilibrium at high temperatures. For too high fhi, the 
lepton number violating scatterings wash out most of the asymmetry after it is generated. 

Scanning over the alloowed ranges for the neutrino mixing angles and taking A = 0.22, 
one can find the regions in the mi — m 2 — m 3 parameter space for which the two conditions 
are satisfied 

• |ei| > 1(T 6 

• 1(T 5 <mi< 1(T 2 

We assume for definiteness in our numerical results that the neutrino Dirac mass matrix is 
identical to the up quark mass matrix[]. There is an undertermined phase in this procedure 
which we can assume is such that the maximal CP-asymmetry is obtained since we are 
determining the points for which it is possible to obtain the required baryon asymmetry. 



4 This parameter is especially important in the super symmetric scenarios there exists a large num- 
ber of scattering diagrams which are lepton number violating, and the Yukawa interactions are 
much more important. 

5 Relaxing this choice, or choosing a different sign for the neutrino mixing angles will change the 
quantitative results, but not the qualitative ones. 
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As an example the allowed regions for Texture 4 is shown in Fig. 1 with the neutrino 
mixings set to s 23 = 0.55, s i2 = 0.07 and si 3 = —0.03. For a particular choice of the 
neutrino mixing angles only a narrow three dimension region is allowed in the full parameter 
space. This region is characterized by larger values of the lightest right-handed neutrino mass 
(10 7 < Mi < 10 8 GeV), and smaller values for the neutrino quark mass matrix parameter 
{rrv D mD)i\- This represents a moderately fine-tuned solution which can be understood from 
the hierarchical structure of the mass matrices as follows. 
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Fig. 1. The allowed solutions for mi and r/1.3 that generate a sufficient baryon 
asymmetry for Texture 4. 
Consider a generic matrix exhibiting the hierarchy given by the textures 

/AX 8 BX % CX A \ 



X 



generic 
D 



B'X 6 DX 4 EX 2 , (14) 
\C'X 4 E'X 2 1 j 

where A, B, B' . . . E' are coefficients of order one. The leading A dependence of the param- 
eters most important to leptogenesis for this generic hierarchy is shown in Table [IT]. Also 
shown are the A dependence of the five RRR textures; one sees that Textures 2, 4 and 5 have 



CP violation e x of the same order as the generic case, while Textures 1 and 3 are further 
suppressed by powers of A (the suppression by A 4 assumes that the texture zeros are exact). 
Furthermore Texture 1 has a enhanced rhi which implies that for this texture the dilution 
factor is very large. 
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TABLE II. The leading A dependence of various parameters important to leptogenesis. Textures 
1 and 3 have leading e\ that is subleading in A, so any solutions are exceedingly fine-tuned. In 
addition, Texture 1 typically has too large a dilution factor because m\ is enhanced. Only Textures 
2, 4 and 5 offer solutions that are not overly fine-tuned. 
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If one assumes the hierarchy suggested by the RRR textures then there must be some 
fine-tuning of the light neutrino masses to get a solution with a sufficient amount of leptoge- 
nesis. That is because the parameter rhi which governs the dilution of the produced lepton 
asymmetry is too large for typical values of the light neutrino masses, related via the seesaw 
as 



*> ~ ^ ■ < 15 > 



The RRR textures fall into the category of the generic texture defined above, where the 
eigenvalues are in the ratio 1 : A" 4 : A -8 . We find that in this case there is some fine-tuning 
required to achieve the required amount of leptogenesis; in Fig. 1 there is only a small pencil- 
like region which produces an adequate lepton asymmetry. For other choices of the neutrino 
mixing angles there is a different linear correlation between the neutrino masses, but the 
same fine-tuning is required. This is because fh\ predicted by the RRR textures is typically 
too large by a factor A 4 and the CP-violation e± is too small by the same factor. If instead a 
more modest hierarchy for the neutrino Dirac masses is assumed, say 1 : A -2 : A -4 , then one 
finds that a typical value of the light neutrino masses allowed by the solar and atmospheric 
neutrino experiments can generated the required level of leptogenesis. This suggests that 
the leptogenesis occurs more naturally in cases with the reduced hierarchy, and the RRR 
textures must be fine-tuned to achieve the required leptogenesis. 

Since leptogenesis can only occur after the end of inflation, the subsequent thermal pro- 
duction of massive gravitinos can occur. The gravitinos interact weakly and the late decays 
of these can modify the observed abundances of light elements or overdose the universe M. 
The right-handed neutrino mass Mi in our numerical solutions is sufficiently light to admit 
a solution to the gravitino problem |TDI . 
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IV. BOLTZMANN EQUATIONS 



The size of the lepton asymmetry that results can be calculated using the full set of 
Boltzmann equations [[TTj]. These have been studied in the scenarios where there are only the 
Standard Model particles |]12|| , but have become available recently in the full supersymmetric 
case as well ||13|| . We consider the supersymmetric case here since the supersymmetric Yukawa 
interactions are sufficient to produce a thermal population of right-handed neutrinos after 
reheating (the nonsupersymmetric model requires the introduction of new interactions [|17| ). 
The Boltzmann equations become quite involved for the supersymmetric case, where it is 
know that the dilution factor can be enhanced over the non-supersymmetric case because of 
the enhanced effect of the Yukawa interactions. 

In principle what is desirable is to scan over all possible values for the masses and mixings 
of the neutrinos that are consistent with the solar and atmospheric neutrino oscillation 
experiments, and to determine the viable parameter choices. We do not do that here for 
three reasons: (1) computational power is exhausted after a few points, as each solution 
of the Boltzmann equation for a parameter choice involves numerically integrating a set of 
differential equations each of which involves a further numerical integration (this integration 
is needed to calculate the reaction density for the two body scatterings which can occur over 
the full kinematic range), (2) the exact equality between the neutrino and up quark Dirac 
matrices is probably only approximate, so our results must be considered qualitative only, and 
(3) there is an unknown phase in the Dirac neutrino mass matrix that controls the amount 
of CP violation in the heavy neutrino decays, so one can only determine an upper bound on 
the amount of lepton asymmetry generated. So we confine ourselves here to demonstrating 
that a particular parameter and texture choice can produce a baryon asymmetry consistent 
with observed result 

Y B = — = (0.6 - 1) x 1(T 10 , (16) 

s 

where ub is the number density of baryons, and s is the entropy density. This quantity con- 
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veniently is insensitive to the dilution that comes about from the expansion of the universe. 
Similar densities Yi can be defined for all number densities n.j. 

Figure 2 shows the evolution of the neutrino densities and the lepton asymmetry as a 
function of the temperature T through z = M\/T for Texture 4 with neutrino masses of 
mi = 1.5 x 1CT 5 eV, m 2 = 3.0 x 1CT 3 eV and m 3 = 4.0 x 10~ 2 eV. This point is one of 
the allowed solutions shown in Fig. 1 that satisfies the requirements for neutrino oscillations 
and for the requirements on e\ and fh\. For these masses the right-handed Majorana mass is 
Mi = 2.9 x 10 7 GeV. Assuming a maximal CP-violating phase, the amount of CP-violation 
from the decays of the lightest Majorana neutrino is e\ = —2.1 x 10 -6 . These masses 
are consistent with the constraints from the solar and atmospheric neutrino experiments in 
Eq. |j. The evolution of the densities proceeds to the right as the temperature of the universe 
decreases. The figure shows the equilibrium density of the lightest Majarona neutrino 
along with the computed density . Nonzero asymmetries of lepton number from fermions 
Yh s and from scalars Yl s develop, and change sign (hence the dip in the figure), and finally 
asymptote to a constants for values of z — M\/T > 3 — 4. Sscattering processes involving 
exchange of supersymmetric particles enforce that Y^ f ~ Yl s - Finally, for completeness, we 
show the scalar neutrino asymmetry for the supersymmetric partner to the lightest Majorana 
neutrino F]_ = Y^o — Y^. This asymmetry also changes sign before eventually vanishing for 
large values of z. The total density of scalar neutrinos Yi + = Y^ c + Y^ is indistinguishable 
from Yn 1 and is omitted from the figure. 
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Fig. 2. The evolution of the fermionic and baryonic lepton asymmetries. The 
asymmetries Yj, / and Yj, s asymptote to a constant value which is recycled into 
a baryon asymmetry that is sufficient to account for experiment if = Yl s = 
-(0.9 - 1.4) x 1(T 10 . 

The baryon asymmetry is related to the lepton asymmetry (in the supersymmetric case) 



via 



Y B = --Y L , 



(17) 



so the observed baryon asymmetry is generated provided that the asymptotic value (for small 
T) of the lepton asymmetry is 



Y Lf = Y Ls = -(0.6 - 0.9) x 10 



-10 



(18) 



So Fig. 2 shows a consistent solution of Texture 4 that explains the baryon asymmetry of 
the universe which is arising after a fine-tuning of the neutrino masses. 
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V. CONCLUSION 



We studied the possibility that the baryon asymmetry of the universe could result from 
lepton number violating decays of heavy Majorana neutrinos. We assumed the Dirac neu- 
trino texture was given by the set of Ramond-Roberts-Ross textures with five zeros (which 
gives rise to correct masses and mixings of the charged fermions). The heavy Majonana 
neutrino mass matrix is obtained by inverting the type-I see-saw formula where the con- 
tributions from the left handed triplet Higgs are neglected, and using the neutrino masses 
and mixings required by the solar and atmospheric neutrino oscillation experiments. The 
lepton asymmetry is produced due to the lepton number violating decay of the lightest right 
handed neutrino. Contrary to naive expectations, the lightest eigenvalue of the heavy Majo- 
rana neutrino mass matrix is in the range 10 5 — 10 7 GeV even though the right handed gauge 
symmetry breaks at M x = 10 16 GeV. This is due to the hierarchy of the Dirac-type neutrino 
texture. We obtained the following results for the feasibility for each texture for generating 
the required baryon asymmetry: (a) A generic neutrino Dirac mass matrix with eigenvalues 
in the ratio 1 : A~ 4 : A~ 8 can produce the observed baryon asymmetry via the baryogenesis 
via leptogenesis scenario in narrow ranges of light neutrino masses. This predicts a strong 
correlation between the light neutrino masses, but only because the masses of the neutri- 
nos must be carefully tuned to achieve the required magnitude of leptogenesis. A neuriino 
Dirac mass matrix eigenvalues in the ratio 1 : A~ 2 : A -4 naturally gives a lepton asymmetry 
of the required level, (b) Textures 2,4,5 generate the amount of leptogenesis expected in 
models with a neutrino Dirac mass hierarchy with eigenvalues in the ratio 1 : A~ 4 : A -8 , (b) 
The position of the texture zeros in Textures 1 and 3 result in a further suppression of the 
generated lepton asymmetry. 

We carried out a detailed analysis of the generated baryon asymmetry by solving the 
Boltzmann equations for a supersymmetric model numerically for Texture 4. This demon- 
strates that the required texture can be compatible with baryogenesis via leptogenesis for 
some specific values of the light neutrino masses consistent with observations in solar and 
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atmospheric neutrino experiments. 
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